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ABSTRACT 


The Solar Electric Propulsion thrust subsystem statistical error model 
was developed as a tool for investigating the effects of thrust subsystem 
parameter uncertainties on navigation accuracy. The model is currently 
being used to evaluate the impact of electric engine parameter uncertainties 
on navigation system performance for a baseline mission to Encke's Comet 
in the 1980s. 

The data given here represent the next generation in statistical error 
modeling for low-thrust applications. Principal improvements include the 
representation of thrust uncertainties and random process modeling in terms 
of random parametric variations in the thrust vector process for a multi- 
engine configuration. 
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I. INTRODUCTION 


The use of continuous thrust for deep space missions presents naviga- 
tion problems that are ordinarily nonexistent in ballistic missions. These 
problems arise from the presence of proportionately large random accelera- 
tions resulting from uncertainties in the thrust vector. Therefore, a princi- 
pal concern in the study of thrust tolerance on navigation accuracy is the 
construction of a realistically adequate model and statistics for the expected 
behavior of the thrust variations. 

In this report, a model for the covariance of random variations in the 
thrust vector is developed in terms of the parameters that affect thrust 
magnitude and pointing. Principal sources of error include the engines, 
thrust vector control and attitude control systems. Section II gives an over- 
view development of the covariance equations and establishes a foundation for 
the technique employed. The structure of the covariance model is motivated 
in brief by summarizing some of the more important aspects of the propaga- 
tion of the covariance of state in the presence of correlated acceleration 
(process) noise. In particular, it is shown that covariance propagation is 
greatly simplified by (1) the decoupling of the mapping from parameter space 
to state space and (2) the use of the Markov assumption in random process 
modeling. Section III summarizes the statistical model and presents a com- 
parison of the results with covariance models used in previous low-thrust 
navigation studies. 

In Sections IV to VII, the details of the calculations required to reach 
these results are presented. In particular, in Section IV, error models are 
developed relating parameter variations to thrust variations. The errors are 
discussed in terms of the physical processes that affect their absolute magni- 
tude and variations as a function of time. The system covariance model is 
derived in Section V. The development considers a general class of error 
sources to be composed of a combination of independent and common errors. 


JPL Technical Memorandum 33-607 


1 



Independent errors, by definition, do not share a common source, whereas 
common or additive errors are derived from a central source. It is shown 
that error sources must be differentiated in this fashion for a multi-engine 
configuration. The distinction is made because the effect of an independent 
source is reduced by averaging, whereas the effect of common sources is 
compounded by addition. The stochastic process for the random process 
variations is given in Section VI. Justification for the use of the Markov 
assumption is discussed in light of expected random behavior of the param- 
eters. In Section VII, the error model parameters are classified in terms 
of independent and common sources. The covariance for the stochastic pro- 
cesses of parameter error sources are subsequently mapped by the error 
models to covariance for component thrust acceleration variations. 

II. DYNAMICS OF STATE COVARIANCE PROPAGATION 


In order to develop an understanding for the structure of the statistical 
model, it is helpful to establish perspective by presenting some of the aspects 
of the dynamics of covariance propagation. The basic considerations are 
simple enough to be understood without recourse to extensive development. 
Reference 1 provides an in-depth analysis of the subject. 

The dynamics of the spatial state vector (y) due to acceleration errors 
is given by a linear perturbation equation of the form 


6y = B6y + 


B = 


_9f 

9y 



> 


where, along the nominal trajectory. 


(1) 



2 
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and where (T,m) represent the nominal thrust and mass respectively. The 
rocket equation 


T 

is used to expand 6(T/m) 



= mv^ (see the Appendix) 

about a nominal point (0), resulting in the relation 



where 

m = propellant mass flow rate, assumed constant 

= exhaust velocity 
A 

Tq = unit vector in the direction of the nominal thrust 

tg = time measured from launch 
rt . 

m(t) = M - m dt, M = launch weight 
For the worst-case condition at encounter, 


mt. 


m. 



0 . 01 


where 


m 


m^ = 


total propellant mass 
spacecraft dry mass 


The quantity 6v e /v eQ is on the order of 6 T/T q so that, to first order, 
the total effect of mass variation and uncertainty can be neglected. A linear 
mapping of thrust acceleration parameter errors to equivalent spacecraft 
thrust acceleration errors is given by 



C.. 

ij 


(N, x Q ) x. 


i = 1,3 

j = 1,8 


( 2 ) 


JPL Technical Memorandum 33-607 


3 



where N is the number of operating thrusters in a multi-engine configuration 

and x„ is the vector of perturbations in the totality of parameters of the 

0 1 
thrust vector process about some nominal value. 

Sufficient accuracy in the statistical description of the thrust process 
is obtained by assuming the thrust parameter errors to be generated by a 
first-order Gauss-Markov random process of the form 

x., = o^x^ + i = 1,8 (3) 

where a. is a constant for the process and is white driving noise charac- 
terized by moments 


Eh-] = 0 


E[n 2 ] = <r 2 

L 1 J Tj. 

1 


The physical justification for the choice of this process will be given in Sec- 
tion VI. 

The covariance matrix for the acceleration noise as a function of thrust 
process parameter noise is the covariance of Eq. (2): 

R (t) = C R (t) C T (4) 

w x 

Typically, an autocorrelation function is 


2 -a. t-T 

R (t. t) = o' e 
x. x. 

l l 


■*The eight components of the parameter vector are listed in Tables 1 and 2, 
where the rotational component due to star tracker errors has been neglected. 
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R X (t) = *1 (t = T ) 

i i 

where cr^ is the process variance. 
x i 

The parameters x are assumed to be statistically independent to first 
order, and the mapping matrix C is such that the component thrust accelera- 
tions are uncoupled; i. e. , 


* = [v x <] 


where the k's and l ' s represent thrust magnitude and thrust pointing sensi- 
tivities respectively. These simplifications are physically justifiable and 
reduce the complexity of computing the covariance for the process noise and 
the consequent propagation in the covariance of state errors. This effect is 
shown by the following calculations. 

Due to the nature of C and the assumed statistics of x, it follows that 
the autocorrelation matrix for thrust acceleration errors must be of the form 


C = 


L k i» * * * . k 5 J [OJ 
[ 0 ] 


*1 0 


o L 


with 


R w.. (t ' T > = 
11 


2 

<r e 
w. 
l 


■a. 

l 


t- T 


i = x,y,z 


where 


R w (t» T ) = 0 for all i f j 

ij 


2 _ V i2 2 . , , 

- 2k. <r . i * 1,6 

x 1 


2 . 22 . 
r w ~ ^ ^ i v. • i" 1,2 

y. z x i 


(5) 
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Through the use of Eqs. (5), random processes for thrust acceleration 
errors can be constructed of the form 


w = Aw + £ 


( 6 ) 


where 


A.. 

ij 


a ( i = j = x, y, z 
i ( 0 otherwise 


and £ is the vector of white driving noise characterized by 



E 






(7) 


The state of the spacecraft is considered to be optimally determined by 
filtering the observables (measurements) in accordance with the minimization 
of some cost function, e. g. , weighted least squares, minimum variance, etc. 
A required input to the algorithm is the propagation in the covariance of state 
due to the corruption by the covariance of acceleration noise w(t). 

The propagation in the covariance matrix of state errors is found by 
the solution of the first-order matrix differential equation derived from 
Eqs. (1), (6), and (7): 

= MZ + ZM T + S 
dt 

Z(0) = z 0 

where 

z = [y, w ] (9 vector) 

Z = E[z z T ] (9 X 9) 
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Parameter matrices S and M are given by 


S 


"(0) 

; (0) 
i 

M = 

B < 
1 

(0) 


~ r 

- - r 


_(°) 

i Q _ 


jo) i 

A 


(9 X 9) 


III. SUMMARY 


A. Strategy 

Previous navigation studies of continuous thrust applications have 
invoked the Markov assumption. However, the resulting autocovariance 
matrix of thrust acceleration errors was assumed to be spherically distrib- 
uted, of the form 


R(t, t) = o-2 e‘°'^ t " T ^[I] (3 X 3) (8) 

2 

where <r is the process variance. The selection of values for the model 

w 2 

parameters (t^, a) were usually made by physical intuition and independent of 
any considerations of thrust parameter variations. 

Although the simplified model given by Eq. (8) is tractable and possesses 
a certain degree of physical justification (Ref. 2), the parameters of the 
model fail to correlate uncertainties in thrust to uncertainties in the dominant 
thrust subsystem parameters and the factors that influence them. To over- 
come this deficiency, knowledge of dominant error sources that contribute to 
thrust magnitude and pointing error was obtained by investigating the composi- 
tion of the component parts of the thrust system and all related subsystems, 
i. e. , attitude and thrust vector control. 

A realistic statistical (covariance) model is therefore achieved by meet- 
ing the following general objective: to specify the parameters of the model 
given by Eqs. (5), i. e. , (oyor)^ z in terms of the totality of parameters 
( x k’ x f) that contribute to errors in the thrust vector. This objective is 
accomplished by the following approach: 

(1) Define a hardware model for the thrust system and all related 
system functions. 
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(2) 


Derive a perturbation model relating changes in thrust in terms 
of engine parameter and related errors. 

(3) Translate expected random behavior of parameter variations into 
statistical variables. 

(4) Map statistics of parameter variations into statistics for thrust 
variations using the thrust perturbation model. 

B. Results 

The behavior of the standard deviations for the normalized components 

of R (t, t), that is, ( r (x, y, z) in percent as a function of time, is shown in 
w ' ' w 

Fig. 1; component autocorrelations are stated in Eq. (15). These data com- 
prise a compilation of thruster performance data presented in Table 1 and 
celestial reference data presented in Table 2. Standard deviation o" Wx is 
shown to be directly proportional to l/\/N and hence follows the power curve. 
Switching points were based on an 18-kW thrust system with six operating 
thrusters. Equal throttling of all operating thrusters was presumed, with a 
maximum throttling constraint of 50% (see Fig. 2b). Standard deviations 
cj- do not follow the lA/N law directly because of the contributions from 

Wy ^ ^ 

cr , where T represents the vector of celestial reference errors mapped 
€ y> z 

to spacecraft body coordinates. 

C. Conclusions 
Conclusions are as follows: 

(la) The standard deviation of the thrust magnitude <r Wx is a root sum 
square of the uncertainty in the totality of parameters that regu- 
late nominal thrust production (see Table 1). The variation in 

c r w with time is attributed to the fact that all independently 
x 

derived error sources vary as n/N on the basis of total thrust 
(ATq) and decrease in proportion to 

<v/N _ _J_ 

N \fN 

on the basis of percentage change in total thrust (AT/Tq). 
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(lb) The maximum value of °' Wx is 3. 5%, which represents a factor 
of 3. 5 increase in the value of (1. 0%) for current process 

noise models in interplanetary flight. Launch and near- 
encounter values of <r Wx approach the value for current models. 

(2) The standard deviations for tr w fluctuate about the 1. 0% level 

y> z 

on the average and represent the uncertainty in thrust pointing. 
These errors are due to the mechanization of the celestial refer- 
ence cr £ and thruster misalignment cr^. Thruster misalignments 
are independently derived, whereas celestial reference errors 
share a common source. Hence celestial reference errors will 
accumulate in proportion to N on the basis of total thrust (AT q ) 
but will be independent of the number of operating thrusters on 
the basis of percentage change in thrust. 

(3) Attitude control limit cycling contributes a negligible amount (less 
than 0. 05%) to thrust pointing uncertainty. This is due primarily 
to the continuous nature of the attitude control; the engine cluster 
(see Fig. 2a) is translated and two of the engines are rotated to 
maintain near-zero position error about the center of mass rela- 
tive to the celestial reference. Limit cycling occurs due to 
backlash in the translator actuator and has a correlation time of 
hours. However, neglecting the contribution from attitude control, 
correlation times for all significant components of the accelera- 
tion process noise is on the order of weeks. 

This result is contrasted with previous process noise models 
which have correlation times for pointing error on the order of 
hours. This assumption reflected the time variation associated 
with the normal period of oscillation of a typical Mariner-type 
deadband attitude control system. 

Consideration of short correlation times is particularly important 
because variations on the order of hours tend to corrupt ground- 
based data types whose normal variations are proportional to the 
Earth's rate of rotation. 

Hence it is concluded that angular acceleration process noise from 
a linear attitude control system should not have a detrimental 
effect on ground-based data types. 
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IV. ERROR MODELS 


A. Definitio n of Fundamental System Functions 

The model for the fundamental system of thrust vector error sources is 
shown in Fig. 2a. The thrust system receives conditioned solar power from 
N power conditioner units. The thruster array is composed of N operating 
thrusters and M spares (a five-thruster array; one spare is shown without 
loss of generality). The switching network couples the power conditioner 
units to the thrusters so that all operating thrusters are connected to separate 
power conditioner units. 

Attitude and thrust vector control (TVC) are achieved simultaneously by 
means of a 2-degree-of-freedom translator mechanism. Control about the 
third axis normal to the plane of translation is achieved by differentially gim- 
balling two of the engines. The TVC mechanism is actuated in discrete steps 
by means of a stepper motor. Attitude control is maintained by acquiring 
celestial references- -traditionally, the Sun and a convenient reference star. 
Sensors used to implement the celestial reference system are two single - 
axis Sun sensors and a star tracker. 

B. Thrust Subsystem 

Before proceeding to the development of the error model for the thrust 
subsystem, a brief description of an electrostatic thruster is given in order 
to promote a qualitative understanding of the physics of electric propulsion. 

A typical thruster is shown in Fig. 3. Thrust is produced by the electro- 
static acceleration of heavy charged particles (mercury ions) through a large 
potential Vg (typically, 3 kV), according to the "rocket equation" (see the 
Appendix, Section I): 


T = mv 

e 

where m is the mass flow rate of mercury (Hg) and v g is the exhaust veloc- 
ity of the ions relative to the spacecraft. 

The ions are formed by the following process, known as electron bom- 
bardment. A source of electrons is provided by introducing mercury into 
the cathode vaporizer C. A baffle prevents high-energy electrons from being 
injected directly into the plasma. The low-energy electrons escape from 
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the cathode cavity and are accelerated toward the anode, which is maintained 
through closed-loop regulation at voltage V^. The electrons spiral outward 
on the magnetic field lines, thus enhancing the ionization efficiency by per- 
mitting a single electron to make many collisions with neutral Hg atoms in 
traversing its path to the anode. 

Mercury atoms are provided by introducing liquid Hg into the main 
vaporizer M. The flow rate of Hg is regulated by sensing the current drawn 
by the screen and accelerator grid power supplies (V , V ) since m is 

S 3, 

directly proportional to the current drawn by the main beam I_ (see the 
Appendix, Section II). The exhaust plume is neutralized by injecting elec- 
trons directly into the exit beam. 

The plasma consists of electrons, singly and doubly ionized mercury 
Hg and Hg , and neutral atoms Hg. The geometry of the electrostatic field 
lines is always such that the particles at tlie plasma interface removed from 
the center of the beamlet diverge as they exit into the exhaust plume (Fig. A-2). 
The resulting divergence causes a thrust degradation per beamlet proportional 
to the cosine of the average value of the divergence angle across the beamlet. 

The nominal thrust from a single thruster is given by the relation (see 
the Appendix, Section II). 


/pj + n/Z r, 2 \ 

T = K \~i' + 2qJ X B /^(cos G cosp)^ 


where 


K 


' m , 


m 0 - mass of atomic mercury (3. 34 X 10~ 25 kg) 


e 

r B 

m 


B 


= electronic charge unit (1. 6 X 10" 19 coulomb) 

= ion current in the exhaust beam = (q + 2q ) m 

1 2 m Q 

- mass flow rate (kg/s) 

= net ion acceleration potential 

= mass fraction of the total flow rate existing as singly and 
doubly charged mercury atoms, respectively 
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cos G = exhaust beam divergence factor, abbreviated Cq 
£ = thrust recovery factor = 1 + e 

(3 = accelerator and/or screen warpage and misalignment angle 

Normalized perturbations in the thrust vector can be written in the 
convenient form 


AT 

T 


_1 

T 


'8T a- -u 
Ax + 

dx 


aT 

8 ? 



( 9 ) 


where T is the magnitude of the thrust and, nominally, T - T x x. The 
appropriate partial derivatives, when substituted into Eq. (9), provide 



where the subscript (0) refers to the nominal value of the parameter. 


2 x, y, z are the vehicle -body- axis coordinate system of unit vectors where 
x passes through the centroid of the exhaust plume and y is aligned with the 
solar panels, as shown in Fig. 2a. 
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The following summary data (Ref. 3) represent current knowledge of 
the uncertainties and contributing factors concerning the parameters given in 
Eqs. (10): 

(1) Vg is the net potential difference experienced by the ions formed 
in the thruster from their point of formation to their point of 
departure from the spacecraft field of influence. This voltage 
will be uncertain to within about 4 V because of varying line drops 
and uncertainties in the thruster plasma potential and the ion 
beam exit potential. An additional voltage uncertainty is caused 
by the regulation of the main beam power supply. For the pres- 
ent units, this is 1%, or 20 V. The combination of these two 
factors gives an uncertainty of about 0. 5% in the thrust and the 
specific impulse. 

(2) Ig is the difference between the currents drawn by the main 

beam and accelerator power supplies and is the main control 

parameter for regulating thrust. Uncertainties in I arise 

B 

primarily from the gain of the control loop, which regulates I , 
and from the drift in the reference that sets Ig. Present regu- 
lation schemes use type 0 controllers. The uncertainty intro- 
duced by the finite gain of the control loop coupled with the 
uncertainty in main vaporizer characteristics is on the order of 
0. 5 to 1%. Also, Ig is set by an analog reference signal against 
which the measured value of Ig is compared. Electronic com- 
ponents used to generate analog signals are subject to thermal 
and time-dependent drifts, which, if uncompensated, can result 
in an error of several percent. With reasonable compensation 
schemes, it is felt that this reference drift can be held to about 

1%. Thus the total uncertainty in I_ is estimated to be on the 

B 

order of 1. 5% rss. 

(3) The electrostatic and mechanical geometries of the accelerating 
structure produce an ion beam composed of many hundreds of 
small, diverging beamlets. The angle of divergence of the 
individual beamlets varies across the exit grid and also varies in 
time as a function of the beam current density. Because of the 
difficulty in measuring an individual beamlet, no precise 
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information is available on the true average beam-divergence loss. 
Faraday probe measurements in the ion beam are generally used 
to estimate the angle of divergence, but the errors in translating 
these measurements into a value of cos 6 are probably large. 

Current estimates of Faraday probe data will not give the diver- 
gence angle to better than ±5 deg. The error or uncertainty that 
this introduces is obviously a function of the angle, which, in 
turn, is dependent on the electrostatic geometry. In general, 
divergence will increase with reduced specific impulse; it could 
vary from about 15 deg at 2 kV screen potential to up to 20 deg at 
1 kV. It is estimated that, at 3000 s, the inherent uncertainty in 
the value of cos 0 will be about + 3% around a base value of 0. 96, 
and that cos 0 will vary with Ig, with the magnitude of this varia- 
tion at present unknown. 

The factor is introduced to account for charge exchange and 
erosion effects. Examination of thruster accelerator grids indi- 
cates that most charge exchange ions originate downstream of 
the accelerator grid. The fast neutrals formed in the process 
then exit with a velocity higher than that of the ions, because 
they have not been decelerated through the full deceleration 
potential. This represents a slight thrust enhancement. A fur- 
ther small thrust enhancement is obtained by the release of 
material from the accelerator grid because of the charge exchange 
ion impact. 

By far the most important factors contributing to subsystem per- 
formance uncertainties are those affecting the mass flow rate. 
Because at present no direct measurement of mass flow rate is 
available, it must be controlled from some a^ priori calibration. 
Present control schemes utilize the relationship between the dis- 
charge power and the mass utilization efficiency (as indicated by 
the ion beam) at constant flow rate to regulate propellant flow. 

This implies that an a priori calibration of (Bj + versus P,pj_j 
(conditioned power delivered to the input terminals of a single 
thruster) and a subsidiary calibration of versus (arc dis- 
charge voltage) are made. In flight, P^^j and are controlled, 
and t\j and are assumed to follow the calibration curves. 
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The difficulty with this scheme is the sensitivity of the calibration 
to a number of thruster parameters, including thruster geometry, 
magnetic field strength and geometry, division of flow between 
main and cathode vaporizers, cathode-keeper potential, total 
extraction voltage, and neutralizer coupling potential. These 
parameters will vary in time as a function of component aging, 
line and load variations, and subsystem random perturbations. 
Using present control schemes, the uncertainty in the initial cali- 
bration is probably on the order of 1%, and the variation in time 
on the order of ±5%. 

(6) The angle (3 represents the achievable alignment accuracy of the 
thrust vector to the nominal thrust direction. This accuracy is 

a function of mechanical tolerances and the thermal load unbal- 
ance on the accelerating grids. No accurate measurements of (3 
are available. However, a careful design should render p less 
than 2 deg (3 cr). Current data indicate that p is time-invariant, 
implying that the grid plates warp to some maximum angles and 
attain a permanent set over the power profile. For the purpose 
of this study, p is considered time-varying with long correlation 
time. 

(7) Thrust vectoring (gimballing) two of the thrusters (see Fig. 2a) 
reduces the net thrust in proportion to the cosine of the gimbal 
angle 6; the maximum gimballing angle is ±10 deg. Because the 
thrusters are used to achieve closed-loop attitude control, no 

a priori prediction of the vectoring loss can be made. However, 

6 is not modeled as a random variable in the thrust error equa- 
tion because it is assumed that the gimbal angles will be cali- 
brated and measured to an accuracy such that the resulting error 
in thrust will be less than 0. 10%. Table 1 summarizes the cur- 
rent best estimates of the various error sources and their 
expected behavior as a function of time in accordance with the 
assumed process noise model. 
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c. 


Celestial Reference Error Model 


Thrust pointing error caused by errors induced in mechanization of the 
celestial reference system can only be estimated in an order-of-magnitude 
sense at this time because sensor mechanization and strategy for acquiring 
the celestial reference have not been established. However, the following 
considerations indicate a probable approach to be taken. 

Ballistic missions in the ecliptic plane require a star tracker with 
aperture-center axis normal to the vehicle roll axis. Consequently this axis 
is directed toward the south ecliptic pole. Since Canopus is the brightest 
star near this location, it is used for roll reference. However, for devia- 
tions from ecliptic flight, the large solar arrays can inhibit the field of view 
because the axis of the solar panels is constrained to be normal to the Sun 
line for the chosen baseline mission; gross rotations of the vehicle about an 
axis tangent to the plane of the orbit thus render Canopus viewing impossible 
during certain portions of the trajectory. The solution of the attitude refer- 
ence problem will possibly involve any one or all of the following considera- 
tions: electrical and/or mechanical gimballing of the star tracker, use of 

multiple reference stars, and Sun- sensor gimballing because of the need for 
solar-panel articulation. A simplified analysis of thrust pointing error 
caused by celestial reference mechanization is presented here in lieu of any 
design data concerning the above technique(s). 

It is assumed that the celestial reference system utilizes the Sun and 
a convenient reference star. The Sun sensors collectively have the equiva- 
lence of a two-axis sensor. Consequently, uncertainties in the Sun- sensor 
output trace a solid angle as shown in Fig. 4. The star tracker is sensitive 
to motion about or out of a plane. Therefore, uncertainties in the star- 
tracker output trace a wedge. 

It is also assumed that (e^, t^) and represent small rotation vectors 
derived from errors in the Sun sensor and the star tracker respectively. In 
terms of percentage change in a single thrust vector normal to T, 

AT _ 1 9T A - 

t ~ t ar Ae 
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where 



( 11 ) 


Pointing error 


can be determined using 


y z w J 

because the drift represents the dominant component of the angular uncer- 
tainty from a statistical point of view. Sensor biases and/or null offsets can 


be estimated in flight quite easily. However, the increase in knowledge of 
the average behavior of the random time-varying component is very slight. 
Nevertheless, the ability to estimate the state of the random process is 
increased as the correlation time becomes large. Sun-sensor errors are 
specified in the coordinate system defined by x', y', c, as shown in Fig. 4. 
The inclusion of star-tracker errors in the model gives rise to thrust cross- 
correlation terms in y and z. However, the effect is probably negligible 
under the given set of assumptions concerning star -tracker drift as shown 
below. 


Table 2 summarizes current estimates of Sun-sensor and star-tracker 
errors and their expected behavior as a function of time. 

D. Thrust Vector Control 


It has been shown that closed-loop operation of the thrust vector control 
system results in a low-frequency, low-amplitude limit cycle in each axis of 
the control system. This oscillation is a function of the electronic compen- 
sation time constants and the amount of backlash in the gimbal and translator 
actuators. 


Laboratory tests have been run to show the effects of a worst-case 
actuator backlash on the system. Table 3 summarizes the results (see 
Ref. 4). Conclusions from the laboratory tests, taken from Ref. 4, are: 

(1) Excessive actuator backlash causes an excessive number of step 
motor pulses in the closed -loop control system. Assuming that 
the total number of lifetime pulses allowed is 10 , it is necessary 
to limit the pulses to one every 3 s on the average. The six-step 
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case shown in Table 3 promotes 20 pulses per cycle, which is 
excessive in terms of lifetime considerations. 

(2) Small-amplitude, low-frequency oscillations are to be expected 
in each axis of the thrust vector control systems. The oscilla- 
tions are so small in amplitude (less than~0. 005 deg peak to 
peak) that no noticeable effect in thrust vector pointing accuracy 
will be observed. 

(3) The mechanization of the actual hardware shows a reduction in 
the total number of step motor pulses by a factor of 3, when com- 
pared to computer simulations of the control system limit cycle 
behavior. 

On the basis of the preceding data, it is concluded that thrust pointing 
error caused by thrust- vector -control mechanization errors is negligible, 
under the stipulation that contributions to thrust pointing error of less than 
0. 0 5% are neglected. This criterion also considers the correlation time with 
equal weight, since variations on the order of hours tend to corrupt ground- 
based data types whose normal variations are proportional to the Earth's rate 
of rotation. 


V. SYSTEM COVARIANCE MODEL 

Variations in the net thrust vector with respect to a set of reference 
body coordinates are considered to occur from two processes: (1) changes in 

the nominal thrust in body coordinates and (2) rotations of the body coordinate 
system. Further classification of the error sources in a statistical sense 
provides data concerning the correlation of random variables. Errors of a 
specific type (i. e. , voltage, current, etc. ) which do not share a common 
source will be independent in a statistical sense, resulting in a root-sum- 
square net contribution. Conversely, an error that shares a common source 
will be perfectly correlated in a statistical sense, resulting in a summation 
of the contributions of the errors individually. However, the collection of all 
thrust parameter variables is assumed to be statistically independent to first 
order. 
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Symbolically, 


W. 

1 


N M 

EE 

k=l j=l 


N L 

A ijk X jk B ijk Y j 

k=l j=l 


for each component of the normalized thrust error vector in the refer- 
ence coordinate system. The summations are carried out over the total num- 
ber of operating thrusters N, the number M of independent errors X., and 
the number L of common errors y . In vector formulation, the thrust error 
is 


N 


N 


W 


= E A i s + E b i 


i= 1 


i= 1 


where, to first order in the errors. 


A - J- i>J- R _ J_ 9T 
A ~T _ 

o ax ay 


( 12 ) 


and where T(t) is the nominal thrust vector per thruster, and Tq is the total 
thrust. Normalized random errors X and y assume the general form 

\ = X 0. + h {t) 

l 


Vi = v 0 + Yj(t) 

l 

where Xq^, Yq. are biases and X^(t), y(t) are time-varying components. It 
is assumed that the magnitude of the error is not a function of thrust and that 
bias and time -varying components are independent. 

It has been assumed that the probable behavior of a single engine is 
representative of the engine cluster as a whole. This assumption provides a 
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great computational savings because only the average response of the 
navigation process to a representative set of uncertainties in the thrust 
process need be investigated. The assumption of statistical independence for 
the elements of the random parameter vectors X and y is entirely justifiable 
for normal operating conditions. 

The autocovariance matrix is the expected value of Eq. (12) at arbitrary 
times tj and t^: 


N 

R w<‘l’ *2» = Z) A i E [VT X i (t 2>] A i T 
i=l 

+ S B i E [v(tj) Y T <t 2 j]^ B 1 T 

i=l i= 1 


(13) 


Define 


where 


A.(t r t 2 ) = E^.ttj), X. T (t 2 )J 
r(tj, t 2 ) = E^y(tj), y (t 2 )J 



E [y i Yjl = o, i ^ j 


Equal throttling of all operating thrusters is assumed, so that the mapping 
matrices A and B remain invariant from thruster to thruster (T Q = NT), 
where 


A(i)’ B (l) = NA » NB 
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Finally, it is assumed that the autocovariance for the X... are identical in 
agreement with a previous assumption, so that 


A. 

1 


R \^ t l’ t 2 ) 


( 0 ) 


( 0 ) 


R \ (W 

M 


for all 


1- ■ ■ N 


Similarly, the autocovariance for the correlated error sources is 


r 


R Y (t,.t 2 ) (0) 

1 • 

(0) R v (t.,t 7 ) 

Y L 1 


Equation (13) takes the form 


R 

w 


(t 


t 2 ) 


n A (i) A ^rV A 


(i) + B d) 


r(ti,t 2 ) 


B 


( 1 ) 


(14) 


The preceding calculations indicate that all independently derived error 
sources (X.) increase in proportion to \JN on the basis of total thrust AT and 
decrease in proportion to <s/N/N = 1 /n/N on the basis of percentage change in 
total thrust AT/Tq. Common errors y add in proportion to N on the basis 
of total thrust and contribute the same amount on the basis of percentage 
change. 

VI. STOCHASTIC PROCESS CONSIDERATIONS 


Choosing a random process to adequately represent the stochastic 
behavior of the thrust and related system parameters is difficult because of 
the lack of any statistical data derived by experimentation. Meaningful experi- 
ments are virtually impossible because of the long correlation times involved. 
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In lieu of this data, some intuitive assumptions about X^(t) must be made. It 
is desirable that y^t) possess the following properties (Ref. 5): 

(1) The process should possess a unimodal probability density func- 
tion. This implies that small values of the noise are expected to 
occur more often than large values. 

(2) The process should be unbiased; i. e. , the statistical average of 
the noise should tend to zero. 

(3) The process should be autocor related in time. This is necessary 
because dominant variations in the process behavior are expected 
to occur at frequencies within the bandwidth defined by the char- 
acteristic frequency of the spacecraft dynamics. 

(4) The process should be stationary. This implies that the variance 
of the noise is expected to remain constant in time. 

A process which fits the preceding description was introduced by 
Ornstein and Uhlenbeck (Ref. 5) as a model for the velocity of a particle 
undergoing a Brownian motion. The statistical properties of the Ornstein- 
Uhlenbeck (O. U. ) process are defined by the following relations; 

(1) The probability density function is unmodal; 

f[x(t)) = a 

V5T x 

where cr is the standard deviation of the process, 
x 

(2) The O. U. process is unbiased; 

E[x(t)] = 0 


(3) The O. U. process is exponentially autocor related in time: 


P ,, , \ 2 “"V*! 


where 1/a is the correlation time of the process and, since 
R (t^ t 2 ) depends only on the time difference (t^ - t^), x(t) is 
stationary. 
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There exists a duality between the continuous Gauss-Markov process 
and the O. U. process because both processes satisfy a Langevin equation of 
the form 

x(t) = a X(t) + u(t) 


where u(t) is Gaussian white noise, i. e. , 

E[u(t)] = 0 

R u (t r t 2 ) ~ Q 6 ( t 2" t ]) 

VII. COVARIANCE MAPPING 

Covariance data for the thrust system and celestial sensors are mapped 
into covariance for percentage change in thrust AT/T by Eq. (14). The 
parameter transformation matrices A and B are 


A 


(1) 


” K 

! <of 



0 ; o 


i 


B d> = 

1 

(0) 

i l p_ 

3X8 

l 

(°) : l 

i __ 


where the elements of the partitions are specified by Eqs. (10) and (11). All 
thrust parameter errors are assumed to be independently derived; thrust- 
pointing error caused by celestial sensor errors must appear as a common 
source to all thrusters. Parameter covariance matrices in normalized form 
are given by the following relations: 


A(t.T) 
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r (t, t) 


e 


t- T 


( 0 ) 


( 0 ) 


y 

0 c r 


- a 


where cr , is the sun sensor variance, and the correlation times of the pro- 

y 

cesses for the thruster parameters are assumed identical without loss of 
generality in accordance with the data of Table 1. Mapping of A and F by 
A and B provides the component autocorrelation functions of percentage vari- 
ation in total thrust for each of the spatial coordinates (x,y, z); 


R (t, t) 
w 

x 




2 2 1 2 1/22 

N l : b c e 4 V B 400 V i 2 . 


40 


1 2~1 

, 000 


, ? -a t-T 

R w <*’ T » * I'p' Xl 


-a t-T 

x 


, -a t-T 

,2 . 2 . e l 

+ cr sin e 


(15) 


R (t,T) 
w 

z 


- o’ ft e 


1 2 -“xM, 2 

+ cr e 
e 

y 


n u P 


An example of thrust error covariance in the spacecraft y axis, caused 
by sun-sensor errors, is shown in Fig. 5. 

The behavior of the standard deviations for the components of R w (t,T) 
as a function of time is shown in Fig. 1. This summary is a compilation of 
Eqs. (15), Table 1, and Fig. 5. Standard deviation cr w is directly propor- 
tional to 1/iv/N and, hence, follows the power curve. Switching points were 
calculated based on an 18-kW thrust system with six operating thrusters. 

Standard deviations cr w do not follow the 1 /\fN law directly because of the 

y » z 

contributions from tr e 

y, z 
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Table 1. Thruster performance summary data 


Parameter 

Nominal 

value 

Calibration 
accuracy 
(% nominal) 

A priori 
process 
standard 
deviation 

Correlation 

time 

Contribution to AT/T 

Bias, % 

Time -varying 
component, % 

X B 

Programmed 

a 

O 

-H 

±1.5% 

Weeks 

0.5 

1.5 

V B 

1 - 2 kV 

±0.5 

±1.0% 

Weeks 

0.25 

0.5 

cos 0 

0.96 

±2.0 

±3.0% 

Weeks 

2.0 

3.0 

^1 

0.8 - 0.85 

±1.0 

±5.0% 

Days -weeks 

0.02 - 0.05 

0.1 - 0.2 

^2 

0.04 - 0.07 

±20.0 

±25.0% 

Days -weeks 

0.5 - 1.25 

0.5 - 1.0 

€ 

0.005 

±30.0 

±30.0% 

Days -weeks 

0.15 

0.15 

P 

0.0 

-- 

±2/3° 

Months 

-- 

1.2 cross axis 



Table 2. Summary of estimated Sun-sensor and star-tracker errors 




Drift 

Sensor 

Bias, deg 

Standard 3- 
deviation, deg 

Correlation 

time 

(tj, € ^ ) Sun sensors 

<0.08 

0.17 

Days - weeks 

(e^) Star tracker 

<0.05 

Negligible 

— 


cl 

Assumes a zero mean uniform distribution between ±0.30 deg. 


Table 3. Effects of actuator backlash 


Actuator 

Backlash, 

steps 

Limit cycle 
period, min 

Limit cycle 
amplitude, deg 

Translator 

1 

30 

0.002 

Gimbal 1 

0.6 

30 

0.003 

Gimbal 2 

6 

30 

0.03 
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percent’ 



Fig. 1. Standard deviations of component normalized 
thrust acceleration process noise 
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Fig. 2a. Fundamental system functions 
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Fig. 2b. Typical power curve and switching strategy 


JPL Technical Memorandum 33-607 



C CATHODE VAPORIZER POWER SUPPLY AND HIGH-VOLTAGE ISOLATOR 
M MAIN VAPORIZER POWER SUPPLY AND HIGH-VOLTAGE ISOLATOR 
V s SCREEN (MAIN BEAM) POWER SUPPLY 
V Q ACCELERATOR POWER SUPPLY 
V d ARC DISCHARGE POWER SUPPLY 

V B = V d + V s 
L = I - I 
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Fig. 3. Electron bombardment engine 
(electrostatic thruster) 
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CELESTIAL REFERENCE: 


/\ 
b = 


l£x$! 


/\ s\ /\ 
a = b X c 


c* = SPACECRAFT/SUN 
DIRECTION 


★ = REFERENCE STAR 
DIRECTION 




Fig. 4. Celestial and vehicle coordinate system relationships 
and celestial sensor error definition 
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percent 



TIME, days 

Fig. 5. Thrust error covariance in the spacecraft y axis 
caused by Sun-sensor errors 
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APPENDIX A 


THRUST EQUATIONS 


I. THE ROCKET EQUATION 

In some time interval At, the impulse imparted to the fluid mass Am 
ejected from the rocket (Fig. A-l) must equal the impulse imparted to the 
rocket. Accordingly, 


L 


Am (v 


0 


v ) - Am v_ 
e 0 


,J = L (m o 


(Fluid Mass) 


Am) v ' - (m Q - Am) 
(Rocket) 


V °_ 


(A-l) 


Let 

m(t) = mg - Am 

Use of the impulse -momentum principle in conjunction with Eq. (A-l) 
provides 


T At = -Am v g = m(t) (v^ - v Q ) 


Let 


At— dt and v^ - v^ — dv 


resulting in 


T 



= m (t) 


dv 

dt 


(A-2) 


II. THRUST EQUATION FOR AN ELECTROSTATIC ENGINE 

The absolute magnitude of the total thrust from an electrostatic engine 
is, from Eq. (A-2), 


T = m v 

e 
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where is the average velocity of the mercury ions at some exit plane in 
the exhaust plume, which is composed of thousands of ion beamlets (Fig. A-2). 
At some time t, let n^ total ions cross the exit plane where n^ and n 2 
singly and doubly charged ions have velocities v ^ and v^ such that 

n l + n 2 = n Q 

The average exhaust velocity is 


v 


e 


n l 

— v i 
n 0 1 



(A -3) 


The energies associated with the singly and doubly charged ions are 


1 2 

7 m 0 V 1 


= eV T 


(A -4) 


E-, = 


I m 0 V 2 


= ZeV B t 


where 

nig = mass of atomic mercury 
e = charge on an electron 
Vg = accelerating potential 
The combination of Eqs. (A-3) and (A-4) gives 


v 

e 


(ti 1 + n/2V| 2 ) 




(A -5) 


where 
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are the mass utilization efficiencies of singly and doubly charged ions and Cg 
is the average value of the cosine of the divergence angle for a single beamlet. 

The total mass flow rate is 

m = mj + rh^ + rrig (A-6) 

where rhg is the mass flow rate of the mercury required for the cathode and 
neutralizer. Define the current in the main beam to be 


. e • , 2e • 

Id = m i + 

B m^ 1 mg 2 


(A-7) 


where 


(A-8) 


Equations (A-6), (A-7), and (A-8) are solved simultaneously to obtain 



m = 


0 

e (ti x 


B 


+ 2q 2 ) 


(A -9) 


The combination of Eqs. (A-5) and (A-9) gives the desired result: 


/2m /r| + n/2-t| 2 > 

T = „ , A - I„ /V„ c, 


e \ Tjj + 2q 2 / J 'B V B 0 
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V EXHAUST VELOCITY OF THE MASS ELEMENT Am RELATIVE TO THE ROCKET 

e 

v n - v ABSOLUTE VELOCITY OF THE EXHAUSTED MASS ELEMENT Am 
0 e 


Fig. A-l. Test rocket for thrust model 



GRID GRID 
THRUST 


Fig. A-2. Single electrostatic thrust beamlet 
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